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Cn Abstract. Searching for transit timing variations in the known transiting exoplanet systems 

T —H J can reveal the presence of other bodies in the system. Here we report such searches for two 

, 1^ ■ transiting exoplanet systems, TrES-1 and WASP-2. Their new transits were observed with the 
4.2m William Herschel Telescope located on La Palma, Spain. In a continuing programme, three 

C , consecutive transits were observed for TrES-1, and one for WASP-2 during September 2007. We 

■ used the Markov Chain Monte Carlo simulations to derive transit times and their uncertainties. 

ri ' The resulting transit times are consistent with the most recent ephemerides and no conclusive 

O ,. proof of additional bodies in either system was found. 
I < 

Q . Keywords. Planetary systems, stars: individual (TrES-1, WASP-2), techniques: photometric 

1. Introduction 

^-H ' Transiting planets provide a wealth of information about exoplanetary systems. Short- 

^ , term variations in the mid-eclipse times of the transits may reveal the presence of moons, 
trojans or other planets (Holman & Murray 2005, Agol et al. 2005, Ford & Holman 2007), 

f^ \ whereas long-term variations could result from orbital precession (Miralda-Escude 2002). 
This provides further constraints on theories of planetary system formation and evolution, 
as well on theories of planetary atmospheres and their composition. 
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O ■ 2. Observations and Data Reduction 



We observed three transits of TrES-1 on UT 2007 September 12, 15, and 18, corre- 
sponding to epochs E = 152, 153, and 154 of the cphcmcris given by Winn et al. (2007): 



H ; Tc{E) = HJD (2453895.84297 ± 0.00018) + (3'!0300737 ± 0'!0000026) x E, (2.1) 

and one transit of WASP-2 on UT 2007 September 13, corresponding to epoch E — 162 
of the cphcmcris given by Charbonneau et al. (2007): 

Tc{E) = HJD (2454008.73205 ± 0.00028) + (2'?152226 ± 0'?000004) x E. (2.2) 

We used AG2, a frame-transfer CCD mounted at the folded Cassegrain focus of the 
William Herschel Telescope of the Isaac Newton Group, La Palma, Spain. AG2 is an ING- 
designed autoguider head with E2V CCD having a field of view (FOV) 3.3 x 3.3 arcmin^ 
and a pixel scale of 0.4 arcsec/pixel. We used a Kitt Peak R filter in order to minimize 
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the effect of color-dependent atmospheric extinction on the differential photometry and 
the effect of limb-darkening on the transit light curve. 

We observed under nearly perfect conditions on the nights of UT September 12, 13, 
and 18. On the night of UT September 15, we observed under scattered clouds in the 
second half of the night. On all nights, autoguiding ensured the positions of all stars on 
the CCD varied by no more than 3 pixels over the course of each night. We strongly 
defocused the telescope in order to minimize the effect of the pixel-to-pixel sensitivity 
variations, and also to enhance the duty cycle. The full- width at half-maximum (FWHM) 
of stellar images was 7 pixels for WASP-2 and ranged from 8 to 13 pixels for TrES-1. 
For TrES-1 for the three nights, we acquired 5, 10 and 10 second exposures, while for 
WASP-2 we used 7 second exposures. 

We used standard IRAl|jJ procedures for the overscan correction, bias subtraction and 
performing the differential photometry. We did not use flat-fields to minimise the impact 
of position-angle dependent scattered light entering the aperture. Different aperture sizes 
were tried in order to find out the one that produced the minimum noise in the out-of- 
transit data. For TrES-1 for the three nights, the optimum aperture radii were 31, 31, and 
34 pixels, while for WASP-2 it was 30 pixels. The sky background was subtracted, using 
an estimate of its brightness determined within an annulus centered on each star. Because 
of the strong defocusing and other field stars, we set an inner radius of the annulus equal 
to the aperture radius and used the width of 10 pixels. The comparison stars used were 
2MASS: J19041058+3638409, J19040934-f 3639195, and J19040792+3640116 for TrES-1, 
and 2MASS: J20304846-H0627355, and J20305168+0628008for WASP-2. For each night, 
differential photometry was obtained by taking the ratio of the signal of the variable to 
the mean of the comparisons, both normalized by a smooth function of time. 

To estimate appropriate error bars for our data, we used a procedure similar to that by 
Narita et al. (2007). We first fitted the light curves with the analytic formulae of Mandel 
& Agol (2002) to find the differences between the data and the best-fitting model. We 
rescaled the error bars to satisfy x^ /Noof = 1-0, where Ndof is the number of the 
measurements in each light curve. For the three nights of TrES-1 we found that the true 
errors are higher by factors of 2.2, 1.6 and 1.9 respectively than the errors including 
only the Poisson noise, and for the one night of WASP-2 by a factor of 2.2. Differential 
photometry with these rescaled error bars and the best-fitting model (solid line) are 
shown in Fig. [T] 

In order to estimate the amount of the time-correlated red noise, we solved the equa- 
tions given e.g. by Narita et al. (2007): 

al^al+al, (2.3) 

crl^^ + ^l, (2.4) 

where cti is the standard deviation of each residual and aN is the standard deviation 
of the average of the successive N points. We selected N to correspond to about 20-50 
minutes depending on the length of the interval with the correlated errors, a^ is the 
white noise, which is uncorrelated and averages down as (1/iV)^/^, Ur is the red noise, 
which is correlated and remains constant for specified N . We adjusted the error bars by 
multiplying [1 -I- N{ar/<JwYY^'^ ^^"^ used these rescaled uncertainties for the subsequent 
fitting procedure. 

f The Image Reduction and Analysis Facility (IRAF) is distributed by the National Optical 
Astronomy Observatories, which are operated by the Association of Universities for Research in 
Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Figure 1. Relative piiotometry for TrES-1 (left) and WASP-2 (right) is shown with the best- 
fitting model (solid line). The residuals and rescaled 1 a error bars are also protted, offset by a 
constant flux for clarity. 

Table 1. Resulting barycentric transit times of TrES-f and WASP-2 systems with their 
uncertainties given by the 68% confidence limits. 



System Epoch Mid-transit time Uncertainty 

TrES-1 152 BJD 2454356.41417 0.00010 

TrES-1 153 BJD 2454359.44431 0.00015 

TiES-1 154 BJD 2454362.47424 0.00020 

WASP-2 162 BJD 2454357.39254 0.00016 



3. The Model 

We used the Markov Chain Monte Carlo simulations with the Metropolis-Hastings al- 
gorithm (Ford 2005) to estimate the statistical uncertainties of the resulting parameters. 
We assumed a quadratic law of the limb darkening. For TrES-1, we fixed the system 
parameters according to Winn et al. (2007) and solved only for the transit times. For 
WASP-2, we fixed the parameters derived by Charbonneau et al. (2007) and fitted the 
limb-darkening coefficients and the transit time. For each night, we created 10 inde- 
pendent chains, with the length typically 100,000 points in each chain. We checked the 
convergence of generated chains using the Gelman & Rubin (1992) R statistic. 



4. Results 

The resulting barycentric transit times can be found in Table [1] Their uncertainties 
are given by the 68% confidence limits. The resulting — C residuals are plotted in 
Fig. [2] for TrES-1, and in Fig. [3] for WASP-2, using the ephemerides dH]), and ((2:21) . 
respectively. In summary, our data provide no conclusive evidence of other bodies in 
either system but the systematic trend of O — C residuals for TrES-1, found by Winn 
et al. (2007), is worthy of further study. For TrES-1, our 3 new transit times plus 1 by 
Narita et al. (2007) extend the interval covered by data for more than 150 cycles. We 
therefore used a weighted linear fit to the transit times as a function of epoch (weights 
inversely proportional to the squares of rnis errors of individual transit times) to derive 
a more accurate ephemeris for the future transit predictions: 

Tc{E) = HJD (2453186.80625 ± 0.00054) + (3'!0300725 ± 0'?0000027) x E. (4.1) 
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Figure 2. O — C residuals of mid-transit times of TrES-1 system. Left. O — C of TrES-1 including 
the data from Cliarbonneau ei at (2005), (their Table 1), Winn et al. (2007), Narita et al. (2007) 
and our three nights from WHT. They were computed using the ephemeris of Winn et al. (2007). 
Right. O — C in detail indicating no significant transit timing variations. 
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Figure 3. O — C residuals of mid-transit times of WASP-2 system show the resulting transit 
time is consistent with the ephemeris given in Charbonneau et al. (2007). 
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